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Abstract 



We derive a semiclassical expression for an energy smoothed autocorrela- 
tion function defined on a group of eigenstates of the Schrodinger equation. 
The system we considered is an energy-conserved Hamiltonian system possess- 
ing time-invariant symmetry. The energy smoothed autocorrelation function 
is expressed as a sum of three terms. The first one is analogous to Berry's 
conjecture, which is a Bessel function of the zeroth order. The second and 
the third terms are trace formulae made from special trajectories. The sec- 
ond term is found to be direction dependent in the case of spacing averaging, 
which agrees qualitatively with previous numerical observations in high-lying 
eigenstates of a chaotic billiard. 
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I. INTRODUCTION 



With the intention of finding a fingerprint of classical chaos in quantum systems, the 
study of eigenstates has attract increasing attentions both theoretically [fl-pOH and exper- 



imentally grin | 

over the last two decades. Among many others, the statistics of wave- 



function amplitude |^(r)| and two-point correlation function (ip(r')ifj(r")) have been widely 
studied. Recently, the study has been extended to the statistics of the Ricaati variable 
of wave function, and a new universality has been found in chaotic and disordered systems 
pOfl . The amplitude distribution of an eigenfunction in a classically chaotic quantum billiard 



system is conjectured to be Gaussian in the semiclassical limit by Berry 0. The numeri- 
cal studies in very high-lying eigenstates support this conjecture [|10],|TT],|15j . In fact, in the 
far-semiclassical limit, even in scarred eigenstates, this conjecture remains true ||15| . 

The correlation function of the wavefunction reveals some local fluctuation properties 
of the eigenstates. Following Berry's conjecture that the wavefunction for a chaotic 
billiard is a superposition of infinite plane waves with random coefficients, the two-point 
correlation function is found to be proportional to Jo{kd) ( J is the Bessel function of the 
zeroth order, k the wave vector and d the distance between the two points). The numerical 
calculations in high-lying chaotic eigenstates reveal some discrepancy with this conclusion 
both in a quantum system whose classical counterpart is the geodesic motion on a compact 
surface of constant negative curvature (as studied by Aurich and Steiner and in the 
Limagon billiard (as studied by Li and Robnik |TT|). Later on, Srednicki and Stiernelof fH| 
investigated this fluctuation from a statistical point of view. 

In a detailed study of the autocorrelation function of an individual eigenstate, Li and 



Robnik [Tl|] found that the two-point correlation function is direction dependent. Only after 
averaging over many directions, do the numerical calculations agree with Berry's prediction. 
To date, a satisfactory explanation for this discrepancy has been lacking. This motivated us 
to study analytically the two-point correlation function of eigenstates in a chaotic quantum 
billiard in terms of periodic orbit theory. 

The billiard system we are considering is a quantum system of one particle moving in a 
closed domain B with area A and boundary dB. The potential vanishes inside the billiard 
and has infinite value otherwise, V^(r) = for all r G B and V(r) = oo for all r G dB. 
The quantum eigenstate of such a billiard satisfies the Schrodinger equation (Helmholtz 
equation), 

= 0, (1) 

with Dirichlet boundary condition, ip(r) = for r G dB, and is normalized, i.e., 
J A drMT)\ 2 = l. 

In order to define a two-point correlation function we first consider a point r in some 
subset of B, say C, with area given by Ac- We will impose two conditions on C. First, all 
points in C are said to be far from the boundary, dB. Second, we fix Ac such that Ac « A. 
Third, we consider a real vector, d, such that for all r G C, r±d/2 lies within the classically 
allowed region but not necessarily in C. We then introduce a function given by 

r(r A. A F) _ £n Mr + d/2)V n (T - d/2) 
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We also consider the position-smoothed object given by 



™ j a p\ Jc^ £„^„(r + d/2X(r-d/2) 

C(C ' d ' AB> = A*S.I«r)|> ' <3) 

Both definitions correspond to a two-point correlation function (the two points are (r = 
±d/2))) of a group of eigenstates. In the first case, it is a function of a point in C, r, 
an energy window, AE, and a real vector, d; in the second case we perform a position 
smoothing in r over the region C. Let N be the number of eigenstates whose eigenenergies 
lie in AE; all N eigenstates must appear in the sum in equations (0) and (|3|). We also define 
an energy E such that 

E - -AE < E n < E + -AE 
u 2 _ n _ u 2 

for all n in the set of eigenstates. 

Thus, from the outset we consider an energy-averaged (energy-smoothed) object as in |7j 
but may, in the end, consider the limit of small AE such that N — 1. We will also mention 
that the correlation function with energy smoothing has also been widely used in the study 
of disordered systems (for a review please see Ref. and the references therein). 

Historically the autocorrelation function in quantum chaos has revolved around the work 
of Berry. Berry's original definition || is similar to equation (H) but considers strictly one 
eigenstate and imposes a technical but important restriction on the area, Ac- Berry considers 
the function given by 

/ c di# n (r)| 2 

with the area, Ac is mentioned in the following. 

From the outset, the quantity in equation (|j) is chosen to be a function of one eigenstate, 
E n , rather than of an energy window containing one or more eigenstates. As for C, it is 
chosen such that its volume tends to zero in the semiclassical limit ||, i.e. lim^o^c — 
but lim h ^ h / y/Ac = 0. As pointed out by Berry |J, the vanishing of Ac ensures the 
smoothed functions to be 'semi classically sharp', whereas \J~Ac~ vanishes more slowly than 
h so that the oscillatory detail on the scale of the de Broglie wavelength is smoothed away. 
This crucial point also allows one to reduce equation (U) to a Bessel function, J (kd), where 
k = V 2mE / h is the de Broglie wave number. Li and Robnik's numerical calculations are 
based on equation (f|) and the condition given above. 

The goal of this paper is to shed light on the above-mentioned direction-dependence 
problem of Berry's autocorrelation function by (4) working with the closely related defini- 
tions, equations @ and @, in the framework of periodic orbit theory. The introduction of 
an energy window, AE, will permit us to do what we want in a clear manner. Moreover, 
by using standard semiclassical methods, we will be able to calculate a first term given by 
Jo{kd) and a second term which carries all directional dependencies of the vector d. 

The paper is organized as follows. In section 2 we give a general semiclassical formula 
for two-point correlation function. In section 3, the effect of energy averaging is discussed. 
The effect of position and direction averaging is discussed in sections. 4 and 5, respectively. 
The paper is concluded by a summary and discussion in section 6. 
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II. SEMICLASSICAL EXPRESSION OF TWO-POINT CORRELATION 

FUNCTION 



To develop equations (0) or (|3|), let us begin by introducing the Green function in the 
energy representation corresponding to equation ([[]). It satisfies the following equation, 

( E + ^ r ,^j G(r"y;E) = 5(r" - r'). (5) 

Using the series expansion for G(r", r'; E) in terms of the eigenstates of the Hamiltonian: 

G(r'\r';E)=j: rn{r " )Mr ' ) (6) 

where e > 0, it is easy to recast equation (|2]) in terms of the Green function, 

E^(r>;(r>(£ -E n ) = ± [G*(r f , r", E) - G(v", v', E)\ (7) 

Take the integral over energy integral [Eq — AE/2, Eq + AE/2] on both sides we have: 

E^(r')C(r") = ^- [(G*(r',v",E) - G(v",v',E))\ (8) 

where (f(E)) represents an energy smoothing or an integral on the energy window, AE, 
containing the N energy levels, 

1 fBo+AE/2 

AE Je -ae/2 

Thus we have: 

J2 |^n(r)| 2 = [<G*(r, r, E) - G(r, r, E))\ = --3m(G(r, r; E)), (9) 



r( j a ci\ 1 (G(t + d/2, r - d/2; E)) - (G*(t - d/2, r + d/2; E)) 

C(r, d; A_E) = — — (10) 

2z ksm(G(r, r; Ej) 

and the position smoothed object, 

rfr+Ajn 1 Ic dr(G(r + d/2 ' r ~ d/2; E)) ~ (G * (r ~ d/2 ' r + d/2; E)) nn 

C(C, d, AE) = / c dr3m<G(r,r;25)> (U) 

In the semiclassical limit, the Green function is decomposed into a sum over two inde- 
pendent contributions: 

G(v", r'; E) = G DI (r", r'; £) + G w (r", r'; E) + O Q) . (12) 
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The first term Gdi on the right represents a contribution from short-range trajectories 
passing directly from r' to r". The second term Djn represents a long-range contribution (i.e. 
they make one or more bounces from the billiard wall) from classical trajectories beginning 
and finishing at these points. 

Now we shall calculate the numerator and denominator in equations ( |I~0D and ( |TTj) sep- 
arately. The denominator, 

Sm(G(r, r; E)) = Qm{G DI (r, r; E)) + %m(G IN (r, r; E)). (13) 

Given, 

dp 



3mG D: (r, r; E) = -vr J 8[E - H(p, r)]- 



(27Tfl) 2 

and that for a free particle moving inside the quantum billiard, H(p, r) = p 2 /2m, we have 

j ii 

Q m (G DI (v,r;E)) = -^ (14) 

The contribution from long-range trajectories (indirect trajectories) is given by the con- 
ventional sum over classical trajectories starting and ending at r 

m , . ,1 /o ( i ~/ „x . ~ 



G '^- r; £) = rn^W ? 1 Ar exp U s(r ' ri £) " 2 )• (15) 



dt"dt' 



t'=i»=0;r'=r"=r 



Here 7 is the count along the trajectory where the determinant (A^) of second derivatives 
of the action in the pre-factor diverges. t'(t") measures a distance perpendicular to the 
trajectory at the initial (final) point and v'(v") is the velocity of the trajectory at the initial 
(final) point. The definition of the coordinates t will become evident shortly. 

Equations fliCf ) and (|Tl"| ) give the general formulae for two-point correlation function. In 
the subsequent sections we will evaluate these two formulae by different smoothing. 



III. THE EFFECT OF ENERGY AVERAGING 

For the rest of the discussion we will write a subscript, d, for all data referring to the 
trajectory, r ± d/2 — > r =p d/2, and for all data referring to the closed trajectory, r — > r. 
Thus, S = S d (r, r; E) and S d = S d (r ± d/2, r =f d/2; E). We have 



$sm(G IN (r,r;E)) 



rn 



E | Ao |V2 cos (|_- (27 + 1) 



(16) 



Using equation (|14D , the denominator (|13D becomes 



%m(G(r,r;E)) 



2k 2 



(17) 
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By the same procedure we can also calculate the numerator, we denote it as In, 
I n = (G(r + d/2, r - d/2; E) -G*{r- d/2, r + d/2; E)) 



I DI + r 



I7V 
JV 



(18) 



where 



'iV 

7-JiV 

'at 



(G D7 (r + d/2, r - d/2; E) - G* DI {r - d/2, r + d/2; E)) 
(G IN (r + d/2, r - d/2; E) - G* IN (r - d/2, r + d/2; E)) 



(19) 



For a quantum billiard the short-range Green function is given by the usual formula, 



G DI {v\v'-E) 



m 
2itf 



H (k\r"-r'\) 



(20) 



where Hq(x) is the zeroth order Hankel function of the first kind and p = |p| = \/2mE, 
k = p/h. Posing r" = r + d/2 and r' = r — d/2 we get 



G D7 (r±d/2,r + d/2,£) 



2ih' 



H Q (kd). 



(21) 



Note that this Green function depends only on the position coordinates via d 
and satisfies the natural boundary condition. The first term in equation fll8f) then becomes, 



m 



and the second term, 



T IN 



2m 



ih^/2-Rh 



1^ = - ij? (J (kd)) 



E|A// 2 (cos(|-^(2 7 + l))) 



cl 



(22) 



(23) 



Thus, 



IN 



m 



(J (kd)) + Vh^=J2\ A d\ 1/2 (cos 
v2vr cl 



Sd _ 7T 

n 4 



(27 + 1))) 



(24) 



We should point out that to obtain the above result, we have used the system's time- 
reversal symmetry whereby the trajectory, r + d/2 — > r — d/2, has the same classical prop- 
erties entering equation (|18"D as the trajectory r — d/2 — > r + d/2. 

Substituting the denominator, equation (|TT|), and numerator, equation (p4"D, together 
into equation (0) we have finally an expression for equation (0) up to the leading term of 

C(r, d; AE) = (J (kd)} + >/»^f £ |^| 1/2 <cos - ^(2 7 + 1))) 

" v^/f £ |Ao| 1/2 (cos (| - ^(2 7 + 1))) (25) 
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The sum in the second term in the above equation is over all trajectories connecting 
r + d/2 and r — d/2. In the third term it is over all closed trajectories starting and finishing 
at r. The energy smoothing has the subsequent affect of giving preference to trajectories 
with a time of motion satisfying, r < H/AE. Thus as we reduce the value of AE we are 
forced to add an increasing number of trajectories in equation ( |25| ) to maintain numerical 
consistency and accuracy Conversely, by choosing an average over many states we are free 
to choose less trajectories in the sum above, an evident advantage from the numerical point 
of view. 



IV. THE EFFECT OF POSITION AVERAGING 

Now we consider the case where numerator and denominator in equation (0) are to be 
integrated over some small region C with area, Ac- This amounts to integrating equation 
(|T7j) and equation (|18D over C. In order to proceed we will first consider the possibility that 
Aq is either fixed (in h) or goes to zero so slowly in the semiclassical limit that the long-range 
Green function has many rapid oscillations on C. 

With the above-mentioned simplification, the long-range Green function can be inte- 
grated via the method of stationary phase. The integral of interest has the general form, 



J dxj dy £ e ^(r±d/2 >rT d/2; B )/ft_ 
cl 



The stationary phase points in C, labeled by r, are solutions of the equation, 

+ (§) - <*> 

(r±d/2,r=pd/2) / (r±d/2,rq=d/2) 




The solutions of equation fl26|) correspond to points r in C for which there exists a trajectory 
leaving a point r' = r±d/2 and arriving at a point r' = f =pd/2 with equal initial and 
final momenta in x and y. This translates into a search for points r in C affiliated with a 
trajectory which has the following two properties: 

• (a) its path includes two parallel chords 

• (b) on one chord there exists a point, r ± d/2, and on the other a point, r =p d/2. 

The first condition assures that the initial and final momenta are equal. The second condition 
assures that a vector, d, can be oriented between the two chords with midpoint given by f . 
Geometrically, we search for a trajectory whose path includes two parallel chords that lie on 
opposite sides (the perpendicular distance is given by dj_/2) of a third chord also running 
parallel to the two outer chords and lying in C (please see Fig.l). It is not important that 
the two chords lie in C, however, by definition r must lie on the midpoint chord which lies in 
this region. It is also clear that the solutions of equation (|2"6] ) are not isolated; in fact they lie 
on a line, the midpoint chord. We speak of a line of stationary phase points, S, and to treat 
the integration in x and y we introduce the usual system of coordinates to describe motion 
local to a trajectory. Moreover, due to time-invariance a solution corresponding to the path 
r + d/2-^r — d/2 will also be a solution corresponding to the path r — d/2 — > r + d/2. 
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There exist two types of trajectories that satisfy the required conditions (see Fig. ([[])). 
Trivially any periodic orbit which has a chord in C parallel to d will satisfy this condition 
and its stationary phase points then lie on that chord (Fig. la). A second possibility comes 
from non-periodic trajectories (Fig. lb). In both cases v" = v' = s. Having found the 
stationary phase points, one can expand the action to second order in both x and y. On 
the other hand, for the purposes of integrating over C it is most convenient to introduce 
the usual Cartesian set of local coordinates, describing motion parallel and perpendicular 
to S. Let us denote this new set as (s, t) with \d(x,y)/d(s,t) \ = 1. The variable t measures 
a displacement perpendicular to S and s measures a displacement along S. For case (1) 
trajectories the point r lies on the trajectory itself but in the second case it lies between two 
chords that are parallel to one another (and to S) and are separated by dj_/2. 

For both cases (1) and (2), one now proceeds to integrate over s and t by expanding the 
action to second order in t (in both initial and final point), for fixed s, and integrating first 
in t. Expansion of the action to second order in t gives, 

S({s + d||/2, t +d±/2), (s - d\\/2, t - d±/2); E) 

= S((s + d\\/2, d ± /2), (s - d ]{ /2, -di/2); E) 
'd 2 S d 2 S 





^ + ^ + W (+0(t) - (27) 

The first term refers to the action, Sj, for going from the point (s + d»/2, d±/2) to the point 
(s — d\\/2, —dj_/2). It is clearly invariant in s, the coordinate along S. Fixing s at one value 
on S, we can define the linear Poincare map, dP((t',p') — > (t",p")). It is a matrix with four 
elements given by: 

St" = mu8t' + mi2<5p' 

5p" = m 21 5t' + m 22 5p. (28) 



We also have the well known relations [E3 

d 2 S mn d 2 S m 22 d 2 S 



dt' 2 ' dt" 2 m 12 ' dt'dt" rri\ 2 



(29) 



which we use to write the second variation as (inn + ^22 ~ 2)/(2mi2). Again in equation 
( P8|) we are considering variations in t at the initial and final points. 
The indirect part of the Green function is given by |p4|, 



G IN (r ,r;E) = E L ,i exp ( - ( ^ + —^-t ] - ij- ] , (30) 




where 



m\ 2 m{ 2 



D(s) = zV, W(s) = "^y 2 - . (31) 



Again the symbol d refers to classical motion beginning at r ± d/2 and finishing at r =p d/2 
whereas the symbol refers to the special case, d = 0. 
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The integration over t gives rise to 

/+oo 
dt exp 
-oo 



'iW{s)^ 



2vrft 



2ft 



m\ x + 171*22 



— — exp (ijsign (m^ + -2)^ (32) 



Here the sign function is given by; 

sign(x) 



1 x > 
-li<0 



(33) 



Due to the invariance of action and the term m u + m 2 2 the subsequent integration in s 
is particularly simple; from it we get the length of S in C divided by the speed along it, |i>|, 
in other words a time which we denote Tj. This gives 



/ dv(G IN (v ± d/2,r t d/2, E)) = -iVhrnV^ 

JA C 



where 



Therefore, 



d J\m'i 1 + m%2-2 



= (exp (^-^); 



a d = 2 7d + 1 - a(m d u + 771*2 - 2) 



/ dr(G IN (r + d/2, r - d/2, E) - G* IN (r - d/2, r + d/2, E)) 

J Ac 



—iVhm 



1 cl 



i fS d cr d . 

1^ / ^(COS — — 7T J ) 

d ylmfi + m^ 2 -2| v ft 4 



The numerator is thus, 
mAc 



In = 



:*2 



ift 



2ft 



( cos ( y ~ If 77 " 



Similarly, we can also obtain the expression for the denominator 



drSm(G(r,r; £)) 



m„4 c 
"1/?" 



1 + 



ft 



E 



1 d 



•Si fTo 

A^ x /| m o i+m 2 -2| NC ° S U 4 71 



Finally the correlation function is given by, 



C(C,d;AE) = (J (kd)) + h— £ 



rpd 

1 cl 




(cos 






- 2 


V ft 






(cos 




\/l m ll + m 22 


-2 


V ft 



!»)> 



-7T 



(34) 



(35) 



(36) 



(37) 



(38) 



The sum in the second term is over all trajectories satisfying the criteria (a) and (b). It 
carries the directional dependency of d via the choice of trajectories. The sum in the third 
term is over all periodic orbits starting and ending at the stationary point in the averaging 
area of Ac, this term does not depend on the direction. All trajectories must be included, 
however, the subsequent energy smoothing will restrict this sum to trajectories whose time 
scale is r « h/AE. Again our final answer cannot depend on AE but this factor will control 
the contribution of each trajectory in the sum. Trajectories whose time-scales greatly exceed 
r will be washed out by the energy smoothing. Any periodic orbit which has a chord parallel 
to d in C will contribute to this term. However, a periodic orbit which scars a particular 
state whose energy is outside the window, AE, may not contribute. 



V. THE EFFECT OF DIRECTION AVERAGING 



In this section, we will consider the sum in equation (BSD over different directions of 
the vector d. To do this one can picture a quantum billiard with vector d centered at 
(x,y) = (0,0). Rotation of d produces a circle of radius d/2 centered at (0,0). Each 
trajectory in equations (|25|) and/or fl3"8|) begins and finishes on this circle. The initial and 
final point are separated by the diameter d. Most generally, the dynamics defines a map P 
which is a function of initial point q = {x,y), momentum p = (p x ,p y ) and time. The linear 
map {DP) is a 4*4 matrix. 



To enumerate trajectories one constructs a symbolic dynamics ||25|| . This construction 
allows one to assign to each trajectory in equations fl25|) and/or (|38|) , a sequence of numbers 
(e.g. a code (i)) and a smooth function Si, the action of that trajectory. Given a trajectory 
labeled by i which starts at 6* and finishes at 9* + 7r, the idea now is to find a new trajectory 
which runs almost the same pattern but begins at 9 = 9* + d9 and finishes at 6' = 9 + n. 
This matching can be performed by using the linear map, DP (with matrix elements given 
by the original trajectory), substituting 5x = —d sin 959, and Sy = d cos 959, and solving for 
initial (5p x ,5p y ). This method for matching can be continued until the linear map has a 
singular point. Clearly Si, the action built for code i, is a function of 9 for all i. 

Let us now denote terms (1),(2) and (3) in equations (p5|) and/or fl38| ) by f , fi and f 2 , 
respectively. We also define an average over 9 given by 



Clearly f and f% are independent of 9 and can be trivially integrated; they remain of order 
1 and \fh respectively. On the other hand, on averaging fx in the semiclassical limit, only 
points 9* where dSi/d9 = will give useful results. After Gaussian integration, any term 
labeled by (i) in fx, will either be of order H if 9* exists in the domain of integration or zero 
otherwise. Thus in the semiclassical limit a sum over many 9 would semiclassically clear 
away corrections of order y/h which are 9 dependent (e.g. the term fx). Moreover, note that 
9x and #2 need not satisfy 9 2 — 9x = 2tt. 

In the paper JTIJ, numerical work on a single high-lying eigenfunction in the Limagon 



billiard was performed. It was shown there that while Berry's conjecture was not satisfied 
for a fixed 9, a direction average over many 9 (with < 9 < 2it) appeared to give a result 
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which exactly corresponded to that of Berry's conjecture. This latter point would intuitively 
suggest that /i exactly cancels f 2 using 9\ = and 9 2 = 2tt, 

^ r f 1 (o)do=-f 2 =-f 2 , 

Z7T JO 

leaving f =< J (kd) > as the only surviving contribution. But according to our arguments, 
fi is of order h while f 2 is of order \fU. Hence these two terms cannot cancel each other. 
However, since fi is washed away by the direction average, the differences from f Q = f Q =< 
J > will come from f 2 alone. 



VI. CONCLUSIONS 

We have calculated a semiclassical expression for an (energy) averaged autocorrelation 
function over a group of eigenstates within the energy window AE up to the leading term 
of \/h. The final expression equations (^) and (|38|) are our main results for correlation 



functions equations (0) and ([|), respectively. In both cases, the autocorrelation functions 
are expressed as a sum of three terms, the first one is the Bessel function of the zeroth 
oder (Berry's term), the second and the third ones are of special trajectory modulations. 
In the first case, i.e. definition equation (0), the second term (in equation (^)) includes 
the contribution from all trajectories connecting the two ends of the vector d = r" — r'. 
The third term comes from the contribution of all closed trajectories starting and ending 
at r. However, in the second case, i.e. have space averaging (equation fl38|) , things become 
quite different. The second term on the right-hand side of equation ([38]) is the contribution 
from all trajectories satisfying the criteria (a) and (b). This term carries the directional 
dependency of d. There are two kinds of trajectories satisfy these two criteria: (1) any 
periodic orbit which has a chord in the averaging domain C parallel to d (see figure 1(a)) 
makes contribution; (2) any open trajectory (see figure 1(b)) having two segments parallel 
to each other and passing through (not necessarily perpendicular to the vector d) the two 
ends of d in the same direction contribute to this term. The averaging over all directions 



might wash away this term as observed numerically by Li and Robnik [LT|. The third term 
on the right-hand side of equation (|38f) includes all periodic orbits starting and ending at 
the stationary point in the averaging domain C. In principle, the number of trajectories to 
be calculated is infinite, however, the energy smoothing introduced from beginning restricts 
this to trajectories and/or periodic orbits whose time scale is about r pa Ti/ AE. This gives 
an advantage in numerical calculation. 
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FIGURES 





FIG. 1. Schematic illustration of two different kinds of trajectories contributing to the second 
term on the right-hand side of equation (38). (a) Periodic orbit which has a chord in C parallel 
to d = r" — r'. (b) Non periodic orbit which has two segments parallel to each other and passing 
through the two ends of d. It is obvious that the contribution depends on the direction of d. 
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